Collider-independent top quark forward-backward asymmetries: 
standard model predictions 
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We compute, for top quark pair production at the Tevatron and the Large Hadron Collider, the 
collider-independent forward-backward asymmetries defined in [l[ in the standard model at next- 
to-leading order in QCD, including also electromagnetic and weak corrections. 
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I. INTRODUCTION 

The charge asymmetry respectively forward- 
backward (FB) asymmetry - in ti production at the Fer- 
milab Tevatron stands out as perhaps the most promi- 
nent anomaly that the data analysis of this collider has 
left. The observable that has been mostly used is the ti 
rest frame asymmetry whose definition is based on the 
difference Ay = y t — yi between the rapidities of the 
top quark and antiquark, which is invariant under boosts 
along the beam direction, 
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FB 



N(Ay > 0) - N(Ay < 0) 
N(Ay > 0)+N(Ay < 0) ' 



(1) 



where N denotes the respective number of ti events. The 
measurements of this asymmetry by the CDF 0, Q an d 
DO [H Collaborations, which are in excess of the standard 
model (SM) expectations have triggered more 

than a hundred theory papers that ex pla in this anomaly 
by new physics (see, for example, |l5H20i |). New physics 
explanations of the anomalous Tevatron asymmetry of- 
ten predict new related effects at the Large Hadron Col- 
lider (LHC) [2l|, including the observation of new parti- 
cles [22]. As yet, none of these effects have been found at 
the LHC. But, of course, this does not rule out the pos- 
sibility that the Tevatron asymmetry results from new 
physics, telling us that, if anything, this new physics is 
perhaps not as simply modeled as by the hitherto existing 
proposals. 

A closer test of the Tevatron excess is provided by the 
measurement of the ti charge asymmetry at the LHC. 
While the Tevatron FB asymmetry in Eq. ([1]) involves the 
rapidity difference Ay, the definition of the LHC charge 
asymmetry used by the CMS and ATLAS experiments 
employs the difference A\y\ = \y t \ — \yj\ between the ab- 
solute values of the top and anti-top rapidities in the 
laboratory (LAB) frame, 



Ac 



N(A\y\ >0)-jV(A|y| <0) 
N{A\y\ >0) + N{A\y\ < 0) 



(2) 



This definition takes advantage of the fact that valence 
quarks q — u,d have a larger average momentum frac- 
tion than antiquarks q. This leads to a boost of the ti 
system along the direction of the incoming quark. There- 
fore, an excess of top quarks in this direction - that is, a 



FB asymmetry in the center-of-mass (CM) frame of the 
initial partons - leads to more t than i quarks for large 
values \y\ of the (anti)top rapidity, while for small values 
of | y | it is the other way around. Current measurements 
of A c by the ATLAS |23] and CMS [H [H] Collabo- 
rations have found agreement with the SM predictions. 
However, these results are per se not incompatible with 
the Tevatron measurements 0, [26l - [28j . since Afb and 
Ac are different observables that result from a differ- 
ent "weighting" of the "intrinsic" asymmetries A u . Ad in 
uu — » ti, dd — > ti, respectively. (Notice that gg — > ti 
does not contribute to Apb and Ac-) 

A direct test of the Tevatron anomaly has been pro- 
posed [1| that consists in the extraction of the asymme- 
tries A u c i from the measurement of the suitably binned 
asymmetries of Eqs. (TT]), at the Tevatron and LHC, 
respectively, and the subsequent comparison of the re- 
spective results. Their numerical values are nearly the 
same at both colliders, up to corrections that are much 
smaller than the experimental precision; thus their de- 
nomination as "collider-independent". The determination 
of the same quantities at the two colliders could shed 
light on the origin of the Tevatron anomalies and settle 
the apparent tension with the LHC measurements. 

The asymmetries A Ut d can be extracted from the FB 
and charge asymmetries in Eqs. (TT]), because they can 
be written, to a good approximation, as 



A FB 
Ac 



A d F d , 

u + AdFdDd ■ 



(3) 



provided we restrict ourselves to a narrow interval in 
the ti invariant mass m t j. This will be shown below. 1 
The factors F q - which differ at Tevatron and the LHC 
and also depend on the CM energy - can be interpreted 
at leading order QCD as the fractions of qq initiated ti 
events. The D q are "dilution" factors that take into ac- 
count that in ti production at the LHC it often hap- 
pens that the initial valence quark has smaller momen- 
tum fraction than the sea antiquark, thus leading to a di- 



1 In fact, in the derivation of Eqs. ||3}, it is the partonic squared 
CM energy s that has to be fixed, which differs in general from 
m t f . But fixing m t j instead of s, which is required in applications 
to data analysis, is good enough for our purpose. 
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lution of the asymmetry generated at the partonic level. 
Both F q and D q depend on the longitudinal velocity of 
the tt pair in the LAB frame, 



/3 



\p!+pI\ 



(4) 



where E, p z are the LAB-frame (anti)top energy and mo- 
mentum along the beam direction, respectively. On the 
other hand, for fixed s, A u ^d are /3-independent. In prac- 
tice, where a finite m t j interval has to be used instead of 
s, A Ut d become mildly /3-dependent. This /3-dependence 
can be weakened by imposing an upper cut on the trans- 
verse momentum pj, of the tt pair. Hence, for a chosen 
interval < m t i < mjj, A u< j, can be extracted from 
a fit to the distributions Afb(/3) an d Ac(j3) measured 
within this m t j interval, as discussed in [lj, using the 
F q and D q factors computed in the SM, e.g. by Monte 
Carlo. 

In this paper we calculate A u ,d in the SM at next- 
to-leading order (NLO) in the QCD coupling, including 
also electromagnetic and weak corrections. To be pre- 
cise, "NLO" refers in this paper to the computation of 
the numerators in Eqs. ([T]), © to order including the 
electroweak corrections of order a 2 s a. In the next sec- 
tion we derive Eqs. (gj) in detail for the SM at NLO. 
This derivation holds also if there are new physics con- 
tributions to the asymmetries ([T]), ©. In addition, we 
discuss the role of gq contributions, which also lead to an 
asymmetry at the LHC, albeit very small in the SM. In 
section [TTTI we present our numerical results. 



II. DERIVATIONS 

The following derivations apply to the computation of 
the numerators of Eqs. (|T|) , (|2]) to NLO in the gauge cou- 
plings (see above). These numerators receive non-zero 
contributions only from terms in the squared matrix ele- 
ments that are asymmetric with respect to the exchange 
of the t and i momenta. As it is well-known, respective 
contributions da A only arise from the matrix elements 
of qq —> tt(g,j), where q — u,d, and of gq(q) — > ttq(q). 
To NLO in the gauge couplings, the charge asymmetric 
terms do qg are infrared-finite, while for qq initiated tt 
production, the soft-gluon divergence that is present in 
d&qq (virtual + soft) cancels against the corresponding 
divergence in J do qq (tig), and likewise for real photon 
radiation. To NLO in the gauge couplings, the numer- 
ators of Eqs. fTJ), © are free of initial-state collinear 
singularities - i.e., no collinear counterterms are required 
to this order. 

The NLO numerators are denoted by A 1 in the fol- 
lowing. For definiteness, we consider the denominators 
of Eqs. ([I}, ([2]) to leading order (LO) and denote them 
by N°. Yet, alternatively, NLO denominators may be 
used, and the derivations are completely analogous to 
the ones presented here. Quantities without subindices 



imply a sum over all partonic sub-processes, whereas a 
subindex, if present, indicates the corresponding subpro- 
cess. For brevity we label with superscripts F, B the ti 
events with Ay ^ 0, respectively, and with superscripts 
>, < the events with A\y\ ^ 0. 

To NLO, the numerator of Eq. flU receives contribu- 
tions from qq — > ti(g,^f), q = u,d. Top-quark pair pro- 
duction by gg fusion is symmetric, and contributions to 
the numerator of Eq. (JlJ by qg, qg processes are com- 
pletely negligible at the Tevatron. Then, the FB asym- 
metry takes the form 



A FB = 



A, 



l.F 



N, 



1,B 



N° 



+ (it -> d) 



A, 



UF - N 1 ' 3 N° 

J * 1111 - 1 * 1111 



N°- N° 
A U F U + A d F d . 



+ (u^d) 



(5) 



Likewise, the LHC charge asymmetry in Eq. @ can be 
written as 



Ac = 



N l,> _ N l,< 


j. N ' 
H 1 


N° 




N i,> _ N i,< 

J Y uu - 1 Y uu 


N°- 

UU 


N°- 

J Y uu 


N° 


+(u -> d) . 





A/i,< 

ly gu 



+ (u -> d) 



N ,- 



N° 



(6) 



The denominators N° in Eqs. (J5J) and © are the (binned) 
LO QCD tt cross sections at the Tevatron and LHC, re- 
spectively. For ease of notation, we use the symbols N°, 
Ngf, etc., both for the Tevatron and the LHC. The SM 
contributions to the numerator of the ch arg e asymmetry 
from gq sub-processes are rather small [l2[ EH and can 
be ignored. 2 (We will explicitly compute them in the 
next section.) Moreover, it will be shown below that, 
provided we restrict ourselves to a narrow s interval and 
to small values of p§, the differences Nq^ 



iv <?<? 



are 



related to the differences — (where "forward" 

and "backward" refer to the initial quark direction) by 



qq qq 1 



N l,F _ N 1 L B 

qq qq 



(7) 



Here D q is a so-called dilution factor, defined, again for 
a narrow interval in s, by 



D„ 



N(x q > Xq) - N(x q < Xq) 



(8) 



9 N(Xq > Xq) + N(x q < Xq) ' 

Using Eq. ([7]) and neglecting the gq contributions, the 



2 If these contributions were measurable, it would be more ade- 
quate to use Ng q for the relative normalization in Eq. ((6j , in- 
stead of N qq . On the other hand, our choice shows more clearly 
the relative size of the gq asymmetries, compared to qq. 
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charge asymmetry (J5]) can be written as 



A c = 



N l,> _ N l,< 

UU UU 


N°- 

UU 


N°- 

UU 


N° 


N 1 '/ - N 1 - 8 

A * uu J Y UU 


N°- 

UU 


N°- 

1 y UU 


N° 



+ (u^ d) 
D u + (u -> d) 



= A' u FiD u + A! d F' d D d . (9) 

Here, we have put primes on A' q and F' q to emphasize 
that these quantities correspond to the LHC, while the 
unprimed quantities refer to the Tevatron. However, as 
it will be shown below, the asymmetries A' u d are, for the 
same narrow interval in s, approximately equal to the 
Tevatron asymmetries A u _ d defined in Eq. ([5|). 

We will first show the equality between A Utd and A' u d . 
For the latter, the forward and backward directions are 
defined with respect to the initial quark momentum di- 
rection. (Of course this is impossible to tell event by 
event.) Then, we derive Eq. ([7]). Our notation is as fol- 



lows. We denote by Xi, i = 1,2 the momentum fractions 
of the initial partons, and f p (xi) is the distribution func- 
tion for parton p in the proton with momentum fraction 
Xi . The dependence of the parton distribution functions 
(PDF) on the factorization scale [ip is not exhibited. 
The (anti)proton 4-momenta at the Tevatron and LHC, 
respectively, are denoted by -Pi, 2, and da denotes the 
differential cross section of a partonic subprocess which 
includes the corresponding phase-space measure e?$ and 
flux factor. 

In the following, we consider binned asymmetries by 
restricting the partonic CM energy s to an interval 
[sminj s max ]. This is accomplished by a factor 

H = 6(s- s min )0(s max - S) (10) 

in the integrals. 

At the Tevatron, the FB asymmetries A U:d are defined 
with respect to the proton direction. Their numerators 
are, in terms of the proton PDF, 
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dx 1 dx 2 f q (x 1 )f q (x 2 )da 1 [q(x 1 Px)q(x 2 P 2 j] [9(Ay) - 9 (—Ay)] % 
+ / dxidx 2 / ff (xi)/ 9 (a? a )<» x [q(x 2 P2)q(xiPi)] [0(Ay) - 9 (-Ay)] H , 



(11) 



where da 1 denotes here the sum of the qq initiated NLO 
differential cross sections for 2-particle and 3-particle fi- 
nal states, and Ay is the difference of the t and i ra- 
pidities in the CM frame of the initial partons. The sec- 
ond integral corresponds to events where the initial anti- 
quark comes from the proton and the quark from the 
anti-proton, and is much smaller than the first integral. 
(It amounts to a "dilution" of order 10~ 3 in the asym- 
metry.) Choosing s max and s m i n in (fTQ|) close enough to 
each other, the factor T-L fixes s within a suitably nar- 
row interval, in which J da 1 , which are functions of s, 
are nearly constant and can then be taken out of the x$ 
integrals. Dropping the argument in da 1 for brevity, 



N. 



1,F 



N, 



l,B 



J d& 1 [9(Ay)-9(-Ay)] 

x / dx 1 dx 2 f q {x 1 )f q {x 2 )U. (12) 



r 



Notice that Ay is independent of X\ and x 2 . The same 
can be done for the LO denominators, 



N ¥ + N lf = / Kb / dX * fMfMn , (13) 



where da qq is the LO differential cross section for qq — > tt, 
so the asymmetries are 



A q = 



J da 1 [9(Ay)-9(-Ay)] 



J dcr qq 

The numerators of the LHC "FB" asymmetries A' u , A' d 
defined in ((9]) are 



N 1,F _ N l lB 

qq qq 



= J dx x dx 2 fqix^f^x^da'iqix^qix^)] [6(Ay) - 6 (-Ay)] H 

dx 1 dx 2 fq(x 1 )f q (x 2 )da 1 [q(x 2 P 2 )q(x 1 P 1 )] [9(-Ay) - 9(Ay)] U , 



(15) 



where the "forward" and "backward" directions are de- fined with respect to the incoming quark direction (note 
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the opposite signs in the arguments of the 9 functions of 
the second integral). By rotational invariance, the sec- 
ond term is equal to the first one, so we can concentrate 
on the former. Taking again the phase-space integrated 
partonic cross sections out of the Xi integrals, we have 



- ( N 1 - 
2 \ ly qq 



Iy qq 



= J da 1 [8(Ay) - e(-Ay)} 

X / dx 1 dx 2 fq{x\)fq{x 2 )%. 



(16) 

The LO denominators = N°k F + N^ B of A' u d are 



1 (n°- 

2 \ ly qq 



qq 



X dxidx 2 f q (xx)fq(x 2 )'H. 



(17) 



Thus the LHC asymmetries A' u d are, given for fixed s, 
by 



A' 



Jda 1 [0(Ay)-0(-Ay)] 
fda° q 



(18) 



They are equal to the ones at the Tevatron, Eq. (fT4|) 



Next we show under which conditions Eq. ([7]) holds. 
The contribution from qq initial states to the numerator 
of the binned LHC charge asymmetry ^ is 



-^ r «f >_ -^95 < = I dxidx 2 f q (xi)fq(x2)da 1 [q(xiPi)q(x2P2)]H£ 

dx x dx 2 f ij (x 1 )f q (x2)da 1 [q(x2P2)q(xiP 1 )]'H£ , 



(19) 



where the asymmetric terms are selected by the factor 

e = 0(\Vt\-\Vi\)-8(\y f \-\Vt\). (20) 

Here yt, yi are the top and anti-top rapidities in the 
laboratory (LAB) frame, respectively. Using rotational 
invariance, the first integral equals the second one, so 
we can concentrate on the former. We now perform a 
rotation-free boost to the tt rest frame. Using that the 



sign of the difference of the t and t rapidities is frame 
invariant, we obtain, with some algebra, that in the limit 
ofp** -> 0, 

£ -)• [9(xi - x 2 ) - 9(x 2 - xi)][6(Ay) - 9 (-Ay)] . (21) 

Inserting Eq. (|2"Tj) into Eq. IT^|) we obtain that for tt 
events with sufficiently small pj, , 



2\ qq qq I 



dxxdx 2 f q (xi)fq(x 2 ) x {9{ Xl - x 2 )6(Ay) - B(x 1 - x 2 )9(-Ay) - 9(x 2 - xx)9(Ay) 

K x 2 - xx)9(-Ay)} x d^H. (22) 



Again, the factor *H fixes s within a suitably narrow in- 
terval, in which the J da 1 are nearly constant and can be 



taken out of the Xi integrals: 



5 



2 \ qq qq I 



da 1 0(Ay) J dx\dx 2 fq{xi)fq(x 2 )9{x 1 - X 2 )'H 
d^ei-Ay) J dx x dx 2 f q (x 1 )f q (x 2 )9(x 1 - x 2 )T-L 
da 1 9(Ay) J dxidx 2 f q (xi)f q (x 2 )0(x 2 -xi)H 
- I da^-Ay) J dx l dx 2 f q {x l )f q {x 2 )6{x 2 -x l )U. 

I 



(23) 



Now let us define 

/ dXldx 2 f q {x 1 )fq(x 2 )da q q6(x 1 



"q 



x 2 )n 



J dx X dx 2 fq{x 1 )fq{x 2 )daqq'H 

J dxidx 2 fg{xi)fg(x 2 )6{x 1 - X 2 )H 

J dx X dx 2 ]q(x X )jq{x 2 )U 

J dxidx 2 fq(x 1 )fq{x 2 )da q q9(x 2 - X^H 

f dXldx 2 fq(xi)fq{x 2 )da q q'H 

J dxxdx 2 fq(xi)fg(x 2 )9(x 2 -Xl)H 

J dx 1 dx 2 f q (xi)fq(x 2 )V. 



(24) 



section for qq — > ttX . In the next section, we use da qq 
in (fMj) . The cancellation of the da q q in the ratios (|2"4~f) 
works, for fixed s, also to NLO - but not to NNLO be- 
cause all the terms in da qq , which is the sum of the contri- 
butions from the tree-level term, virtual corrections, soft 
and hard gluon radiation, and the collinear counterterm, 
are convoluted with the same product of PDF. 



Clearly, d r q + d™ = 1. The integrals with respect to Xi 
in Eq. (|2"3")1 can be written in terms of d r „ and , resulting 



Here da qq denotes the (LO or NLO) differential cross in 



N±<) = J da 1 9(Ay)d r q J dx 1 dx 2 f q (x 1 )f q (x 2 )n - J <B l 6(-Ay) d r q J dx 1 dx 2 f q (x 1 )f q (x 2 )'H 

-Jda^iAy) d™ J dx 1 dx 2 f q (x 1 )f q (x 2 )H + J da 1 6 (-Ay) d™ J dx^f^f^x^H .(25) 



Rearranging terms, we have 

l^f-iV^) = K-<) J dx Y dx 2 f q (x Y )U(*2) 

xda l \6{Ay)-6{-Ay)]U. (26) 
Comparing with Eq. p^|) we obtain Eq. ([7]), i.e., 

N lt ~ N lt = K ~ O ( N ¥ ~ N qf) > (27) 

where 

d r q - d% = D q (28) 

are the dilution factors introduced in Eqs. ([7]), (JSJ). 

We recall that this derivation holds for fixed s and 
sufficiently small p". On the other hand, the formula (|27|) 
and the resulting formula ^ holds for arbitrary values 
< (3 < 1 of the longitudinal velocity of the ti system. 
In practice, the requirement of fixed s must be replaced 



by choosing a reasonably narrow bin in the tt invariant 
mass m t (, i.e., 

H ^ 0(mu - m L )6( mu - m ti ) . (29) 

This will be done in the numerical computations of the 
next section. Then the intrinsic asymmetries A\ will 
become /?- dependent; that is, the formulae Eqs. (fT4")) and 
(UHJ) do no longer apply - for the computations one has 
to use instead the definitions of A ud given in Eqs. ([5]) 
and ©, respectively. But we will show that, for a given 
m t i bin, this /3-dependence is rather mild in the SM to 
NLO. More importantly, as the results below will signify, 
A' u and A' d remain equal to A u and A^, respectively, to 
a good approximation - even if no upper cut is imposed 
on p". In addition we will show by numerical compu- 
tation that neglecting the gq contributions to the LHC 
charge asymmetry Ac in the formula ([5]) is indeed justi- 
fied, given the level of precision one aims at in applying 
Eq.© to future data analysis. 



6 



III. NUMERICAL RESULTS 

Our numerical calculations are based on the code de- 
scribed in 0, [l4[ • We compute the binned asymmetries 
of Eqs. ([T]) and @, and the asymmetries A Ut d (Tevatron) 
and A' u d (LHC, 7 and 8 TeV) for a sequence of intervals 
to_l < m t t < mjj . Within a specified interval [mL,mjj], 
the asymmetries Afb and Ac are computed, for ft bins 
of width Aft = 0.2 for the Tevatron (i.e., < ft < 0.2, 
etc.) and Aft = 0.1 for the LHC (i.e., < ft < 0.1, etc.). 

In the numerators of Eqs. ([T]) and and of A</) u( j, we 
take into account the 0(cP s ) QCD and the 0(a 2 s a) elec- 
troweak corrections. For definiteness, we evaluate the 
denominators of all asymmetries considered in this pa- 
per with LO QCD matrix elements, which is in the spirit 
of a consistent fixed-order perturbative expansion of ra- 
tios like Eqs. ©, ©. The fractions Fp and the dilution 
factors D q are computed for the Tevatron and the LHC 
(7 and 8 TeV) using LO QCD matrix elements both in 
the numerators and denominators. We evaluate both the 
numerators and denominators of the binned asymmetries 
and of F£ ' and D q with NLO parton distribution func- 
tions. 

As emphasized above, the analysis for A[/ d could also 
be done by replacing, on the left- and right-hand sides 
of Eqs. ([5]) and ((6]), the global normalization factors N° 
(Tevatron and LHC) by the respective NLO factors N 1 . 

We use m t = 173.1 GeV (on-shell mass), the QED 
coupling a(mz) = 0.008, and the weak mixing angle 
sin 2 9 W = 0.23. We use the CTEQ6.6M PDF and 
the respective value of a s (mz) provided by this set. We 
put hr — /if = M> an d numerical results are given for 
/i = rn t /2,m t , and 2m t . These scale choices are purely 
conventional. In Ref. [l| the asymmetries A ( J d were ob- 
tained for a benchmark new physics model using a two- 
parameter fit to the Afb (ft) and Ac(ft) distributions, 
mimicking the procedure that has to be eventually per- 
formed with real data. That can be done, for the Teva- 
tron and the LHC, by minimizing 



E 



[dA FB {fti)Y 
[Ac[fti) - A' u FMfti) - A' d F'MftAf 



[dA c (fti)Y 



(30) 



with respect to A ( u ] and A d . Here i labels the different 
ft bins and dApsifti), dAc(fti) are the statistical uncer- 
tainties of the binned asymmetries. Unfortunately, this 
procedure requires extremely high Monte Carlo statistics 
in order to have the two-parameter fit converging to the 
true values. Especially at the LHC, the /3-binned A' u d 
are obtained from the ratio of a tiny asymmetry Ac(fti) 
over a small F q D q (fti) factor. Therefore, in order to save 
computing time, we calculate the asymmetries with a 
one-parameter fit, considering uu and dd contributions 



separately. The values of A ( u } d presented in Tables Q] - Mil 
below are obtained from a one-dimensional least squares 
parameter fit, 



A n 



A' 



E 



A%{fti) 
[dA%{fti)i 



i [dA^m] 2 




[dAfB^)]- 



,(3i; 



where the superscripts of ApBi Ac indicate that we 
restrict the calculation to these specific sub-processes, 
eventually including gq contributions as well. We will 
demonstrate below the consistency of both methods by 
showing that the values of A u ' d calculated using either 
([3H)l or (|3"Tj) agree very well within the expected experi- 
mental uncertainties. The one-parameter fit is more pre- 
cise. 

We estimate the statistical uncertainty of the /3-binned 
asymmetries by taking an integrated luminosity of 20 
ftr 1 for the Tevatron and 10 (30) fir 1 for the LHC with 7 
(8) TeV. This corresponds to an eventual combination of 
results from both experiments at the Tevatron and LHC, 
respectively. A selection efficiency of ~ 25% is taken for 
the semileptonic ti decay channels, similar to that found 
in the experimental analyses [2^, |24|. This results in a 
combined efficiency factor of about 7%. The results for 
the Tevatron are collected in Table U and for the LHC in 
Tables HTland lllll without and with an upper cut < 30 
GeV. 

TABLE I: Asymmetries A u ,d at the Tevatron, without cut 
(top) and for p" < 30 GeV (bottom). 







M = 


m t 


A* = 


2mt 


M = 


m t /2 


mtt 


[GeV] 


A u 


A d 


A u 


A d 


A u 


Ad 


< 400 


0.058 


0.039 


0.054 


0.036 


0.061 


0.044 


400 


-450 


0.096 


0.066 


0.091 


0.060 


0.102 


0.073 


450 


- 500 


0.123 


0.086 


0.116 


0.079 


0.131 


0.095 


500 


- 550 


0.145 


0.102 


0.137 0.092 


0.154 


0.113 


550 


- 600 


0.164 


0.115 


0.156 


0.106 


0.176 


0.128 








M = 


m t 


A* = 


2mt 


M = 


mt/2 


m t t 


[GeV] 


A u 


A d 


An 


A d 




A d 



< 400 0.069 0.046 0.065 0.042 0.075 0.051 
400 - 450 0.117 0.078 0.110 0.071 0.126 0.087 
450 - 500 0.150 0.101 0.141 0.092 0.161 0.113 
500 - 550 0.178 0.120 0.167 0.109 0.191 0.135 
550 - 600 0.201 0.137 0.190 0.125 0.217 0.153 



These tables show that the SM values of A u ^ and A' u d 
computed for the Tevatron and the LHC, respectively, 
are in quite good agreement, which is remarkable given 
the difference of roughly one order of magnitude between 
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TABLE II: Asymmetries A' u d at the LHC (7 TeV), without 
Pt cut (top) and for p" < 30 GeV (bottom). 



ji = m t fx= 2m t fJ, = mt/2 
m tf [GeV] A' u A' d A' u A' d A' u A' d 
< 400 0.055 0.038 0.052 0.035 0.059 0.042 
400 - 450 0.089 0.060 0.084 0.055 0.096 0.066 
450 - 500 0.112 0.077 0.106 0.070 0.120 0.085 
500 - 550 0.128 0.083 0.120 0.076 0.136 0.092 
550 - 600 0.142 0.093 0.134 0.085 0.151 0.101 
600 - 650 0.155 0.103 0.146 0.093 0.165 0.113 
650 - 700 0.164 0.110 0.156 0.102 0.177 0.122 
700 - 750 0.176 0.119 0.165 0.104 0.185 0.129 
750 - 800 0.182 0.118 0.170 0.107 0.195 0.131 



TABLE III: Asymmetries A' u d at the LHC (8 TeV), without 
Pt cut (top) and for < 30 GeV (bottom). 



m t i 


[GeV] 


M = 


m t 




2m t 


fj, = mt/2 






A u 


A' d 


i ' 


A' d 


A' A' 


< 400 


0.055 


0.038 


0.052 


0.035 


0.059 0.042 


400 


-450 


0.088 


0.058 


0.083 


0.053 


0.094 0.065 


450 


- 500 


0.111 


0.075 


0.105 


0.070 


0.118 0.084 


500 


- 550 


0.126 


0.081 


0.119 


0.075 


0.135 0.091 


550 


- 600 


0.139 


0.089 


0.131 


0.081 


0.147 0.097 


600 


- 650 


0.153 


0.102 


0.144 0.094 0.163 0.113 


650 


- 700 


0.161 


0.107 0.153 


0.099 


0.172 0.120 


700 


- 750 


0.172 


0.115 


0.162 


0.106 


0.185 0.129 


750 


-800 


0.177 


0.115 


0.167 0.106 


0.190 0.131 







M = 


m t 


M = 


2m t 


H = mt/2 


m t t 


[GeV] 


4' 


A' d 


4' 




A' A' 

sl-u ^d 


< 400 


0.071 


0.054 


0.068 


0.047 0.077 0.059 


400 


-450 


0.115 


0.078 


0.108 


0.071 


0.124 0.087 


450 


-500 


0.149 


0.103 


0.140 


0.093 


0.160 0.114 


500 


-550 


0.170 


0.110 


0.159 


0.100 


0.183 0.124 


550 


-600 


0.193 


0.128 


0.180 


0.116 


0.209 0.144 


600 


-650 


0.211 


0.143 


0.197 0.127 0.227 0.159 


650 


- 700 


0.229 


0.153 


0.215 


0.140 


0.247 0.169 


700 


- 750 


0.245 


0.165 


0.228 


0.148 


0.261 0.184 


750 


-800 


0.262 


0.177 0.244 


0.158 


0.282 0.194 



[i = rat A* — 2mt fi = rat 1 2 

mtf [GeV] K A' d A' u A' d A' u A' d 

< 400 0.073 0.040 0.068 0.050 0.078 0.061 

400 - 450 0.115 0.074 0.108 0.070 0.123 0.086 

450 - 500 0.147 0.101 0.139 0.092 0.159 0.113 

500 - 550 0.170 0.118 0.159 0.102 0.183 0.124 

550 - 600 0.191 0.128 0.179 0.113 0.205 0.141 

600 - 650 0.209 0.148 0.197 0.126 0.227 0.157 

650 - 700 0.229 0.157 0.215 0.139 0.246 0.168 

700 - 750 0.242 0.164 0.225 0.147 0.261 0.186 

750 - 800 0.258 0.174 0.243 0.155 0.277 0.193 



the predictions for the inclusive asymmetries Afb and 
Ac- For illustration, in Fig. [T] the intrinsic asymmetries 
in Tables [T MlIII are displayed, for /i = m t , as functions 
of m t t- For the Tevatron, only bins with m t i < 600 
GeV are included. These plots show the increase of these 
asymmetries with increasing m t i and, furthermore, that 
the upper cut on reduces the slight difference be- 
tween the Tevatron and LHC asymmetries, making them 
nearly equal. We remark that the differences exhibited 
in the left plot are irrelevant from an experimental point 
of view, as we shall see below. We also point out that the 
differences between the Tevatron and LHC results origi- 
nate from Monte Carlo statistics to some extent, as it can 
be noticed from the fact that the smooth increase of the 
asymmetries with m t i is modulated by small fluctuations. 

As we have emphasized before, at the LHC the contri- 
butions from gq processes to A' u d are quite small. This 
can be shown, for example, by using in the second of 
Eqs. 431]) the SM values oi A qq+9q fa) instead of Agfa). 
The resulting intrinsic asymmetries are presented in Ta- 
ble [TV] for 7 TeV, without cut. Comparison with the 
corresponding numbers of Table [H] shows that the dif- 
ferences are negligible. The differences are even further 
reduced if an upper cut on is applied. 

Next, we check the equivalence of the two- and one- 



TABLE IV: Asymmetries A' u>d at the LHC (7 TeV), without 
Pt cut, including gq contributions. 



/i = rat = 2mt fJ, = rat/ 2 

mtt [GeV] A' u A' d A' u A' d A' u A' d 

< 400 0.056 0.038 0.053 0.035 0.060 0.042 

400 - 450 0.091 0.060 0.086 0.055 0.097 0.067 

450 - 500 0.115 0.078 0.108 0.071 0.123 0.087 

500 - 550 0.132 0.085 0.124 0.078 0.141 0.095 

550 - 600 0.147 0.097 0.139 0.088 0.157 0.106 

600 - 650 0.161 0.108 0.152 0.097 0.172 0.118 

650 - 700 0.172 0.115 0.163 0.107 0.186 0.129 

700 - 750 0.184 0.125 0.172 0.110 0.195 0.136 

750 - 800 0.191 0.125 0.178 0.113 0.205 0.139 



dimensional fits in determining the intrinsic asymmetries. 
In Fig. [5] these asymmetries are plotted for the first four 
m t i bins. The dots represent the values calculated with 
one-parameter fits (these are the numbers given in Ta- 
bles Hl lIIII) , and the ellipses are the two-dimensional 68% 
confidence level regions from the 2-dimensional fit, where 
the centre is the best-fit value giving the minimum Xmin 
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FIG. 1: Summary of asymmetries for the Tevatron and the LHC (7 and 8 TeV), for fi = mt, without cut (left) and for 
Pr < 30 GeV (right). The data correspond to the first two columns of Tables [TI 4IIII 



and the border corresponds to x 2 = Xmin + 2.3. From 
these plots, it is also clear that the slight differences be- 
tween the Tevatron and LHC asymmetries are irrelevant 
as to the anticipated experimental uncertainty. This jus- 
tifies the ansatz of extracting the same quantities from 
two different sets of data. 

So far we have determined the intrinsic asymmetries 
■^ud k v simulating the proposed fitting procedure with 
SM data: we have computed, for various m t j bins, the 
binned asymmetries Afb(P%), Ac (Pi), the fractions Fp 
and the dilution factors D q in the SM and performed the 
fits using Eqs. (|31|) and (f3T))) . This leads, by definition, to 
constant A { u , for each m t i bin. It remains to show that 
this is an acceptable procedure - i.e., that the /3-binned 
intrinsic asymmetries are only mildly /3-dependent within 
the m t t bins chosen above, as it was claimed. This is 
shown in Fig. [3] for the first m t i bin. (For the other m t i 
bins the behavior is quite similar.) This variation can be 
compared, for example, with an increase in F' u by a fac- 
tor of 3 and D u (Dd) by factors of 20 (40), between the 
bins < j3 < 0.1 and 0.9 < j3 < 1. These results corrobo- 
rate the assumption of constant A[},, especially when an 

upper cut on pj, is used. The LHC results shown in the 
plot in the right panel of Figj3] exhibit some statistical 
fluctuations, which have some effect on the resulting fit 
values of of A' j . 



IV. CONCLUSIONS 

The formulae ^ allow to extract the intrinsic forward- 
backward asymmetries A u% d and A' u d from the Teva- 
tron ti forward-backward asymmetry Afb and the LHC 
charge asymmetry Ac, respectively, if measured in suit- 
ably chosen bins of the ti invariant mass m t i and longitu- 
dinal velocity j3 of the ti system. We have shown under 
which conditions Eqs. © hold in the SM to NLO in the 
gauge couplings. Our derivations apply of course also 
to possible new physics contributions to these asymme- 



tries. In particular, we have shown within the SM that 
the intrinsic asymmetries are indeed collider-independent 
and, furthermore, only mildly /^-dependent for suitably 
narrow m t i bins, especially if an upper cut on the Pj, of 
the ti samples is applied. This corroborates the proposal 
of [l[ to use ^ with constant A u and Ad for performing 
a two-parameter fit to the respective Tevatron and LHC 
data. 



In order to apply Eqs. © to data analysis, one has 
to compute the fractions F q and factors D q in the SM, 3 
either to LO or NLO QCD for a specific PDF set and 
for chosen values of the renormalization and factoriza- 
tion scales. The outcome of the fits depend, of course, 
on these choices. Needless to say, this would be extremely 
valuable information. It has previously been shown, on a 
model-independent basis [l[ as well as for specific new 
physics models (26l - [28| . that current measurements of 
A fb at the Tevatron and Ac at the LHC are compatible. 
Therefore, the comparison of the measured values A B * p d 
with their SM predictions A s " d would reveal whether or 
not there is agreement with the SM - and if a deviation 
would be found, it would show whether it is located in 
A u or Ad, or in both. This would be a big step forward 
in pinning down the origin of the new physics contribu- 
tion^), if there are any. The experimental determination 
of A u and Ad will certainly be a challenge, as it will in- 
volve in general a 3-dimensional unfolding of the data 
with respect to m t t, j3, and the (anti)top rapidity, but it 
is certainly worth the effort, for the aim of resolving this 
puzzle. 



3 The present knowledge about the (differential) hadronic tt pro- 
duction cross section implies that possible new physics contribu- 
tions to these functions can be neglected. 
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FIG. 2: Comparison between asymmetries obtained with two- and one-dimensional fits, for fi = mt, without p T cut. The mass 
range in the lower left corner indicates the m t j bin (in GeV). 
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FIG. 3: Dependence of A Uid on j3 in the SM, for m tf < 400 GeV at the Tevatron (left) and the LHC with 7 TeV (right). Black 
and gray triangles correspond to no p" cut and cut < 30 GeV, respectively. 
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